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Abstract 

It is shown that hadronic matter formed at high temperatures, according to the 
prescription of the statistical bootstrap principle, develops a critical point at nonzero 
baryon chemical potential. The location of the critical point in the phase diagram, 
however, depends on the detailed knowledge of the partition function of the decon- 
fined phase, near the critical line. In a simplified version of the quark-gluon partition 
function, the resulting location of the critical point is compared with the solutions of 
other approaches and in particular with the results of lattice QCD. The proximity of 
our solution to the freeze-out area in heavy-ion experiments is also discussed. 
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1. Introduction 

Critical QCD at nonzero baryon density and high temperature is a theory of fundamental 
significance associated with the properties of primordial matter (quark matter) both in the 
early universe and in the multiparticle environment, created in heavy-ion collisions. The¬ 
oretical arguments, based on the universal features of chiral phase transition, suggest with 
confidence the existence of a critical point in QCD which is the end point of a quark-hadron 
critical line of first order and belongs to the Ising universality class [1]. The location of the 
critical point in the phase diagram needs a more detailed treatment, based on first principles 
and it is only recently that a first attempt to locate the critical point, within the framework 
of lattice QCD, has appeared [2], It is of interest to note that the knowledge of the location 
of the critical point in the phase diagram is of primary importance for heavy-ion physics 
since one expects to observe spectacular critical fluctuations only in those collisions of nuclei 
which freeze-out close to the critical point [3]. 

In this Letter we present a different, complementary approach in the search for the 
existence and location of the critical point, based on the statistical bootstrap principle [4-6]. 
We notice that in the vicinity of the QCD critical point (in the crossover regime) the hadronic 
degrees of freedom are thermodynamically indistinguishable from the quark matter degrees 
of freedom, suggesting that the QCD critical point must also be present in the solution 
associated with the statistical bootstrap equations. These equations describe the dynamics 
of interacting hadronic matter at high temperatures, incorporating in the partition function 
the observed hadronic spectrum [4-6]. 

In the context of the Statistical Bootstrap the interaction among hadrons is simulated 
by the presence of clusters. The bootstrap hypothesis adopts that every cluster is composed 
by a number of other clusters described by the same mass spectrum f(m 2 , 6, s ,...), where 
b stands for baryon number, s for strangeness and the dots for any additional quantum 
number that can characterise a cluster, like charge q, absolute strangeness (number of strange 
plus antistrange quarks) |s|, etc. Then fdm 2 equals to the number of clusters with mass 
m and quantum numbers equal to the set b, s ,.... Every cluster is also characterised by 
a kinematical term B(m 2 ), apart from its mass spectrum which is of dynamical origin. 
The conservation of four-momentum, as well as, conservation of quantum numbers can be 
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accommodated by a suitable bootstrap equation [7] 
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The pattern of clusters composed of clusters continues until the known hadrons, consid¬ 


ered as the “input” particles, are reached. The bootstrap equation after performing a series 
of Laplace transformations acquires the form 1 . 

<e(0, {A}) = 2 G(P, {A}) - exp[G(/9, {A})] + 1 . ( 2 ) 


In the last equation G(/3, {A}) is the Laplace transform of the mass spectrum with the 
accompanying kinematical term 
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and (p(P, {A}) the Laplace transform of the input term 2 
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The masses correspond to the masses of the input particles, which in this paper will 

be all the known hadrons with masses up to 2400 MeVs, the gb qs \s\ are degeneracy factors 
due to spin and the A’s are the fugacities of the input particles. 


The bootstrap equation [7] defines the boundaries of the hadronic phase through the 
relation 


_ 7 >(/MA}) < In4 — 1 . (5) 

1 We shall use the abbreviation {A} for the set of fugacities corresponding to the quantum numbers 
(b,q,s, M), i.e. (X b , X q , X S ,7 S )- 

2 The Statistical Bootstrap Model is formulated in the Boltzmann approximation which is acceptable 
throughout the phase diagram except for low temperatures and high baryon densities, where the effect of 
the Bose-Einstein and Fermi-Dirac statistics becomes important [5]. For example for the temperature and 
baryon density which correspond to SPS the error due to the neglection of the correct statistics does not 
exceeds 3% [7]. 
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The same equation leads to the asymptotic mass spectrum of interacting hadrons 

f(m 2 , {A}) "^-3° 2C({X})m~ a ~ 1 exp [m/3* ({A})] , (6) 

where /3 = T l and j3* corresponds to the maximum inverse temperature [7]. The kinematical 
term B and the dynamical term t can be redefined so that their product remains unchanged. 
The different choices adopted lead to different values of the exponent a in eq. (6) and to 
different physical behaviour of the system through different partition functions. In this Letter 
we shall fix a = 4. This choice leads to normal behaviour of the system near the boundaries 
of the hadronic phase since it does not allow the energy density to become infinite, even for 
pointlike particles [6]. Consequently, for this choice, the maximum hadronic temperature 
is not a limiting temperature, since it can be reached at finite energy density and so it 
is consistent with a phase transition [8]. It also guarantees that the bootstrap singularity 
associated with a quark-hadron phase transition is reached in the thermodynamic limit [9]. 

2. The existence of a critical point 

For the particular choice of a = 4, ip assumes the form 

V>(/MA}) ^ 2 J(3B ^( Aa ({ A » + VIA })~ 1 )Y.g™ m l i K 1 (Pm a , i ) , (7) 

where “a” represents a particular hadronic family and “i” the hadrons in the family with 
different masses. The partition function of the pointlike interacting hadrons for this choice 
of a reads [7]: 

OO OO OO OO 
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b '=—oo q'=~ oo s'=—o o |s'|=0 

4 RV r°° 

—=- / x 3 G(x, {\})dx = Vf(j3, {A}) , (8) 

where the only left parameter B is the energy density of the vacuum (bag constant). In¬ 
cluding corrections due to the finite size of hadrons (Van der Waals volume corrections) one 
takes into account the repulsive part of the interaction [10]. With this improvement, the 
partition function in the grand canonical pressure ensemble, given by the Laplace transform 
of Z(V. /3, {A}), is written [10] 

>r(f./3,{A})= e-zbl+V {A}) , (9) 
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where / is the point-like partition function (8) divided by the volume and the Laplace variable 
£ gives a measure of the external forces acting on the system. In what follows we consider 
hadronic systems free of external forces and therefore we fix the variable £ at the value £ = 0 
[9]. From eqs. (8) and (9) one obtains the energy and baryon number density (e, uj,) and the 
pressure P of the thermodynamic system with the help of the following equations: 

-l 
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The onset of a phase transition can be traced by studying the pressure-volume isotherm curve 
(Fig. 1). It is obtained with the help of eqs. (10) and it reveals a region of instability where 
pressure and volume decrease together. This is a signal of a first-order phase transition and 
is due to the formation of bigger and bigger clusters as the hadronic matter approaches its 
boundaries. This instability is removed by a Maxwell construction as it is shown in Fig. 2 
and discussed in more detail in what follows. As expected, no such behaviour is exhibited if 
the interaction incorporated in the statistical bootstrap hypothesis (SB) is neglected. In fact, 
in an ideal hadron gas (IHG), even with Van der Waals volume corrections, the pressure is 
everywhere a decreasing function of volume, keeping the system far from a phase transition. 
The results are displayed in Fig. 1, where uq is the normal nuclear density = 0.14 fm~ 3 . 
For the ideal hadron gas we have used in place of /(/?, {A}) the function: 

]nZ IHG (V,/3,{\}) 1 


V 


2ii 2 /3 
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( 11 ) 


where g ai are degeneracy factors due to spin and isospin. In both cases (SB, IHG) we have 
imposed the constraints: < S >= 0 (zero strangeness) and < b >= 2 < q > (isospin 
symmetric system, i.e. the net number of u and d quarks are equal). 

The existence of a maximal pressure in the isotherm curves (Fig. 1), obtained in the 
bootstrap solution (eqs. (10)), is the starting point in the search for the QCD critical point. 
In fact the end point of the critical line of a first order transition is reached when two isotherm 
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curves corresponding to the two phases and the same temperature, intersect at the point of 
maximal pressure, on the hadronic side. The point of intersection, in this case, coincides 
with the critical point which in general corresponds to nonzero baryon chemical potential. 
In order to quantify the above description we have to consider the partition function of the 
quark-gluon system (QGP) and study the pressure-volume isotherm curves. For this purpose 
we adopt a minimal description of non interacting quarks and gluons within a bag supported 
by an external pressure B. More specifically, the free energy is written [ll] 3 


(5 1 In Zqgp{V ,, /?, X u , Xd) — 


N S N C V v 


Ip 2 + m? e/V^+^A" 1 + 1 


- BV . 
45 


where the index i runs over all quarks and antiquarks and N S1 N c refer to spin and colour 
degrees of freedom. The current masses are taken rn u ps 5.6 MeV, ~ 10 MeV, m s ~ 
200 MeV and the fugacities satisfy the constraints: X u = A^ 1 , A= Aj 1 , A s = A^ 1 = 1 
(strangeness is set to zero). The first two terms in eq. (12) give the contribution to the free 
energy of the noninteracting quarks and gluons respectively and the third term gives the 
contribution of the vacuum (bag constant). 

In order to study the P — V isotherms of the quark-gluon system, we extract from the 
free energy (12) the baryon number density in quark matter, as follows 


ui q \^X u ,X d ) 


N S N C ^ , Ar f°° p 2 dp 

2tt 2 ^ 1 1 Jo e^ + Xi ' 


(13) 


The index i includes only u, u and d, d quarks whereas N t = ±1 for quarks ( u , d) and 
antiquarks (h, d) correspondingly. For a given value of 5, the point where the hadron and 
quark phases meet can be traced on the P — V isotherms for different temperatures. One 
has to require that at the same baryon density = vj 1 ' 1 (eqs. (10) and (13)) the pressures, 
for the two phases, coincide. As can be seen in Fig. 2, for a low temperature, the quark 
matter isotherm meets the corresponding hadron isotherm at a point where the pressure is 
an increasing function of volume. This region is associated with the onset of a first-order 
3 The set of quark flavour fugacities {X u , Xd, X s ,j s ) can be used equivalently in the place of the set 
(Ab, Xq, As, y s ) [7]. 
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phase transition and a Maxwell construction is needed in order to remove the instability 
[6,9]. As the temperature increases the quark matter isotherm meets the hadronic isotherm 
at a point which corresponds to a higher density and pressure. With this set of intersection 
points, the critical line of a first order quark-hadron phase transition is approximately built 
up, in the phase diagram. At a certain temperature, the two phases meet at the point of 
maximal hadronic pressure (Fig. 2). The Maxwell construction is no longer needed since 
we have reached the end point of the critical line, in other words, we have reached a critical 
point of the strongly interacting matter. For higher temperatures we enter the crossover 
regime where quark and hadron phases are no longer distinguishable. 

In Fig. 2 the limitations of our approach are also illustrated. In fact, in the absence 
of a final theory, based on QCD, which could provide us with a unified treatment of the 
equation of state in both phases, the matching of the isotherms, shown in Fig. 2, introduces 
artificial discontinuities in the derivatives for T > T C} in a region where singularities are not 
allowed (crossover regime). Moreover, the scaling behaviour expected at the critical point 
and expressed through the isotherm critical exponent (5), is not recovered at the point of 
intersection for T = T c (Fig. 2). Despite these drawbacks, our treatment concerning the 
existence and location of the critical point may still remain valid if in the exact theory of 
the equation of state, not only the unwanted discontinuities are washed out but, at the same 
time, the qualitative features of the hadronic isotherm, predicted by Statistical Bootstrap 
(Fig. 1), remain unchanged. 

3. The location of the critical point 

On the basis of the above mechanism and in order to locate the critical point in the phase 
diagram ^ — T, one has to determine the parameters (/?, X u , A d , X s , X' u , X' d ) by solving the set 
of equations: 



(14a) 

P < '‘ ) (ftA„,A i ,A»)=P<«>(/3,A'„,Ai) 

(14b) 

dpW(0,\ u , A„,A,) 
dX u 

(14c) 

( S({3, X u , Xd, X s ) ) h — 0 

(14d) 
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(14e) 

(14f) 


( b(f3, X u , A d, A s ) ) h — 2 ( Q(/3, A„, A^, A s ) ) h — 0 

( b ( P ,\' u ,\' d )) q -2 ( Q ( PX , K)) q = 0 

Equations (14d)-(14f) account for isospin symmetry and zero strangeness. 

For a given value of the bag constant B we have investigated the solution of eqs. (14) 
in order to build up the critical line (excluding eq. 14c) and locate the endpoint (including 
eq. 14c) in the phase diagram. We have projected out the solution onto the plane /p, — T 
and our Erst remark is that for a wide range of values of the bag constant B ( B C 4 < 282 
MeV) the critical baryon chemical potential fi c is not zero (Fig. 3a). One may therefore 
establish, on the basis of the statistical bootstrap hypothesis, the existence of a critical point 
in hadronic matter at high temperatures and nonzero baryon chemical potential. A typical 
solution is obtained if we impose the condition /^(hadronic matter) = /ij,(quark matter) at 
the critical point. The solution obtained in this case corresponds to the value B 1 ^ = 250 
MeV and it is illustrated (SB) in Fig. 3b. The location of the endpoint (critical point) in this 
solution is fixed by the critical values T c « 171 MeV and ~ 385 MeV whereas the critical 
temperature at zero chemical potential (at the end of the crossover) is T 0 ~ 188 MeV. In 
the same figure, the solution of lattice QCD with rif = 2 + 1 (LQCD), corresponding to the 
critical values T c = 160 ± 3.5 MeV, fi c = 725 ± 35 MeV and To = 172 MeV is also shown 
[2] for comparison. The discrepancy of these two solutions is due (a) to the fact that the 
lattice QCD result for the critical point is based on rather unphysical input values of the 
quark masses m u , rrid and (b) to the uncertainties in the partition function (12) coming from 
neglected interactions in the quark-gluon system. One expects that, when both, the quark 
masses get reduced and the quark-gluon pressure for nonzero chemical potential becomes 
available on the lattice [12], the SB and LQCD critical lines in Fig. 3b will converge to a 
unique solution, supported both by QCD and Statistical Bootstrap. Finally, it is of interest 
to note that within a different frame of approximate theories (Nambu-Jona-Lasinio (NJL) 
model or a random matrix (RM) approach) the QCD critical point has been located in the 
region T c rs 100 MeV, ~ 600 — 700 MeV of the phase diagram [1], in a good distance from 
the SB and LQCD solutions (Fig. 3b). 

In Fig. 3b the freeze-out points associated with recent heavy-ion experiments are also 
shown on the plane — T [13]. One observes that the SB solution brings the critical 
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point close to the freeze-out area of the SPS experiments [13]. However, before drawing any 
conclusion concerning the phenomenology of the SB critical point, an improved lattice QCD 
calculation of the pressure of the quark-gluon system with small quark masses is needed [12]. 

In conclusion, we have shown that the hadronic matter develops a critical point at high 
temperatures and nonzero baryon chemical potential, in accordance with the requirements 
of the QCD phase diagram. Our approach was based on the Hagedorn Statistical Bootstrap 
Principle (SB) as it was elaborated, in connection with quark-hadron phase transition, in 
references [6,9]. With the present degree of approximation for the quark-gluon partition 
function, the location of the SB critical point, in the phase diagram, was found close to the 
freeze-out points of the SPS experiments with nuclei (T c « 171 MeV, /i c ~ 385 MeV). 
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Figure Captions 

Fig. 1 The pressure-volume isotherms for hadronic matter described by Statistical Boot¬ 
strap (SB) and for the ideal hadron gas (IHG). 
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Fig. 2 The intersection points of the pressure-volume isotherms of the two phases are shown, 
for different temperatures, and the Maxwell construction associated with the phase transition 
is illustrated. 

Fig. 3 (a) The SB solutions for the critical line and the corresponding critical point, for 
different choices of the bag constant, (b) A typical SB solution is shown in the phase 
diagram (//;, — T) together with the LQCD solution and the freeze-out points of heavy-ion 
experiments. Also the NJL-RM solution for the location of the critical point is shown. 
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